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The semiclassical solution of quantum Dirac constraints in generic constrained 
system is obtained by directly calculating in the one-loop approximation the 
gauge field path integral with relativistic gauge fixing procedure. The gauge in- 
dependence property of this path integral is analyzed by the method of Ward 
identities with a special emphasis on boundary conditions for gauge fields. The 
calculations are based on the known reduction algorithms for functional determi- 
nants extended to gauge theories. The mechanism of transition from relativistic 
gauge conditions to unitary gauges, participating in the construction of this solu- 
tion, is explicitly revealed. Implications of this result in problems with spacetime 
boundaries, quantum gravity and cosmology are briefly discussed. 

1. Introduction 

It is well known that the path integral is an efficient tool for solving the Schrodinger 
equation. Numerous results in modern physics in this or that way are related to 
the path integral method. The virtue of this method is that it is applicable to any 
evolutionary (parabolic) equation and for the Schrodinger equation gives its unitary 
evolution operator. At the same time there exists a class of problems in which the 
fundamental dynamical equations are not of a manifestly evolution type. A wide 
class of such problems is represented by Dirac quantization of constrained dynamical 
systems |], ||. In such systems the Schrodinger equation is supplemented by quantum 
Dirac constraints on quantum states, so that the problem amounts to solving the 
whole system of equations, not all of them being of the evolutionary type. Moreover, 
in the parametrized systems with a vanishing Hamiltonian there is no independent 
Schrodinger equation and their quantum dynamics is encoded in the Dirac constraints 
along with their gauge invariance properties. Applications of the path integral method 
in this context are much less known and generally look as follows. 
Consider dynamical systems with the canonical action 
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in configuration space of canonical coordinates and momenta (q,p) = {q l ,Pi) and La- 
grange multiplyers iV = iV 7 '. The variation of iV 74 leads to nondynamical equations - 
the constraints 

T M ( 9 ,p)=0. (1.2) 

The constraint functions on phase space T^(q,p) belong to the first class when they 
satisfy the Poisson bracket algebra 

{T ll ,T v } = U^T a (1.3) 

with some structure functions U® u = U^ u (q,p). This algebra indicates that the theory 
possesses a local gauge invariance under the action of canonical transformations of 
(q, p) generated by constraints themselves and by certain transformations of Lagrange 
multiplyers [||, |J. 

Dirac quantization of the theory ( |1 . 1| ) consists in promoting initial phase-space 
variables and constraint functions to the operator level (q,p,T^) — > (q,p, T M ) and se- 
lecting the physical states | ^Sf) in the representation space of (q,p, T M ) by the equation 
T^l^f) = HI, H, |, 0. Operators (q,p) are subject to canonical commutation re- 
lations [q k ,pi\ = iTi5 k and the quantum constraints T M as operator functions of (q,p) 
should satisfy the correspondence principle with classical c-number constraints and be 
subject to the commutator algebra 

fa,t] = %hufrr x . (i.4) 

with certain operator structure functions standing to the left of operator con- 
straints. This algebra generalizes ( |1.3| ) to the quantum level and serves as integrabil- 
ity conditions for quantum constraints. In the coordinate representation, (q\^f) = 
Pk — hd/idq k , the latter become the equations on the physical wave function 

fp(q,hd/idq)V(q) = (1.5) 

with the differential operators of quantum Dirac constraints T^q, Tid/idq). 

Note that without loss of generality we did not include in (|1.1|) the Hamiltonian 
H(q,p) nonvanishing on constraint equations ( |1.2| ). By extending the phase space of the 
theory with extra canonical pair (q°,po), q° = t, subject to the constraint p +H(q,p) = 



one can always reduce the system to the case of action fll.lj ). At the quantum level 
this extra Dirac constraint plays the role of the Schrodinger equation for the wave 
function of the theory with parametrized time, *f?(t, q) = *f?{q°, q), 



j^ + H(q,hd/idq) 



*(<Z°,<Z) = 0. (1.6) 



In this sense the dynamical content of any theory can be encoded in the quantum Dirac 
constraints of the above type. 
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The path integral in this context arises as a special solution of quantum constraints 
(|1.5|) in the form of the two-point kernel K(q, q') - the analogue of the two-point 
evolution operator for the Schrodinger equation || |[ 

f,(q,hd/idq)K(q,q') = 0. (1.7) 

Similarly to the theory of the Schrodinger equation this is a path integral over con- 
figuration space variables (q(t), N(t)) in (space)time domain £_ < t < t + with the 
boundary conditions related to the arguments of this kernel q(t + ) = q, q{tJ) = q'. 
However, in view of the non-evolutionary nature of equations (|1.7p the path integral 
construction here is much less straightforward than in the Schrodinger case. 

Mainly the motivation for such applications comes from quantum gravity theory 
H |9|, [TUL [II], [12[]. In this theory ( [TTJ ) is a canonical form of the Einstein action with 
<t = 9ab(x-) - 3-metric coefficients and ~ ((— g 00 ) _1 / 2 (x), g 0a (-x)) - lapse and shift 
functions in (3+l)-foliation of spacetime |13|], ( |1.2| ) are the gravitational Hamiltonian 
and momentum constraints and ( |1.5| ) represents the system of Wheeler-DeWitt equa- 
tions on the cosmological wavefunction *&(q) = [ <7 a &(x) ] ||. The construction of the 
path integral here is persuing two main goals. The first goal is a two-point solution for 
the Wheeler-DeWitt equations (Q B 1, H 0, which in view of parametrized nature 



of time encodes the dynamical information. The second goal consists in specifying the 
distinguished cosmological quantum state - the model for initial conditions in quantum 



cosmology of the early universe [10]. The first formulation of the path integral for the 
two-point solution of Wheeler-DeWitt equations belongs to H.Leutwyler [|J. It was 
only qualitatively correct because at that time the structure of gauge fixing procedure 
and the role of ghost fields in the path integral have not yet been understood. The 
path integral with unitary gauge fixing procedure and exhaustive set of boundary con- 
ditions was later proposed in [[J. Then this canonical path integral was converted to 
the spacetime covariant form of the functional integral over Lagrangian variables in 



relativistic gauges |TT, [12 



Of course, these results encorporate a well-known statement on equivalence of the 
canonical and covariant quantizations pioneered in H. In contrast to this, the works 



H |TD , IT], [12] were focused on the nontrivial boundary conditions in spacetime. Correct 
treatment of these boundary conditions leads to the proof that this path integral solves 
the quantum Dirac constraints ( |1.7| ). However, this proof given in |9|, [10], [TJ|, |II| has 



a formal nonperturbative nature and does not even allow one to fix the operators of 
quantum constraints. One can only infer from this proof that these operators satisfy 
the correspondence principle with their classical counterparts and have quantum cor- 
rections which in a rather uncontrollable way depend on the calculational method for a 



path integral |14j| . Thus, no check of the solution to quantum constraints was thus far 
given by direct calculations of the path integral. The goal of this paper is to perform 
such a check. This will be done for generic systems subject to first class constraints of 
the above type in the one-loop approximation of semiclassical expansion. 
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The nature of this check is certainly a comparison of two results: one obtained by 
solving the equations ( [L.7|) and another by calculating the path integral. In the one-loop 
approximation for the Schrodinger equation such a comparison is based on the Pauli- 



Van Vleck-Morette formula |15j and the reduction method for functional determinants 
||16|| . The preexponential factor in the subleading semiclassical order is given by the 
Van Vleck determinant - the solution of the corresponding continuity equation. On 
the other hand, it is given by the functional determinant of the inverse propagator of 
the theory - the contribution resulting from the one-loop (gaussian) approximation for 
the path integral. The equality of these two expressions follows from the reduction 



method for functional determinants considered in [161. The semiclassical solution of 



quantum Dirac constraints in terms of the modified Van Vleck determinant was built 



in the author's paper [|H], Here we reproduce this solution from the Faddeev-Popov 



path integral of ref . [Ill 



The paper is organized as follows. In Sect. 2 we recapitulate the semiclassical solu- 
tion of quantum Dirac constraints of refs. [IT], ITS]] with a particular emphasis on unitary 



gauges participating in its construction. Sect. 3 serves as a link between the canonical 
formalism of the constrained system and its Lagrangian version with local gauge in- 
variances presented in (space)time covariant form. The corresponding Faddeev-Popov 
path integral in relativistic gauges constructed as a solution of quantum Dirac con- 



straints in |TT[ is presented in Sect. 4 along with its semiclassical expansion up to the 
one-loop order. Sect. 5 deals with the mechanism of its gauge independence based on 
Ward identities and, in particular, gauge independence of its local measure cancelling 
the strongest power divergences. In Sect. 6 we calculate the contribution of the gauge 
field functional determinant by deriving the algorithm of its reduction to the Van Vleck 
elements of the semiclassical solution of Sect. 2. Similar calculations are performed for 
the ghost field determinant in Sect. 7. They accomplish the proof of the main result - 
path integral derivation of the solution to quantum Dirac constraints. In concluding 
section we discuss implications of this result regarding the aspects of gauge invariance 
in quantum cosmology and Euclidean quantum gravity. 

The final remark, which is in order here, concerns our notations. Throughout the 
paper we use condensed DeWitt notations || which allow one in a manageable form to 
handle both the quantum mechanical and field theoretical problems within the language 
of the above type. We imply that the range of indices of q l and 

i = l,...,n, /i = l,...,m, (1.8) 

in field models formally extends to infinite dimensionalities of phase space n and space 
of gauge transformations m, and these canonical condensed indices together with dis- 
crete tensor labels carry also continuous labels of spatial coordinates x (or certain 
discrete numbers when the fields are expanded in the basis of some countable set of 
spatial harmonics). Contraction of these indices will include spatial integration (or 
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the corresponding infinite summation)!]. Starting from Sect. 3 we shall also need the 
covariant condensed indices including the time label and allowing one to represent 
the covariant (space)time operations of integration and differentiation in the form of a 
simple contraction. 



2. Canonical solution of quantum constraints: uni- 
tary gauges 



The solution of quantum constraints ( |1.7| ) found in |T7], 



K(q,q')=P(q,q')e 



(2.1) 



contains the Hamilton- Jacobi function S(q, q') which determines the tree- level approx- 
imation and the preexponential factor P(q, q') accumulating loop corrections. The 
system of Hamilton- Jacobi equations for S(q, q') 



T,(q,dS/dq) = 



(2.2) 



enforces the quantum constraints in the leading order of fi-expansion. For their operator 
realization proposed in []T7], [TJ, Q the subleading (one-loop) order yields the following 
quasi-continuity equations for the preexponential factor P(q, q') (see also [19|] in the 
gravitational context) 



8 

— (V { P 2 ) 

dq^ ' 

BT 



K x p\ 



p = dS I dq 



(2.3) 
(2.4) 



A particular solution of eqs.( [2.2|) - ( |2.3|) corresponds to the choice of the principal 
Hamilton function for S(q,q') - the classical action (LI) calculated at the extremal 
of equations of motion g = g(t \q,q') joining the configuration space points q and q 1 . 
Together with eq. fl2.2p this function also satisfies the Hamilton- Jacobi equations with 
respect to its second argument 



TJq',-dS/dq') = 0. 



(2.5) 



Similarly to the WKB theory of non-gauge quantum systems, this function gives 
rise to a special one-loop preexponential factor |] which is a generalization of 

1 This formal approach does not certainly incorporate a rigorous handling of ultraviolet infini- 
ties and possible quantum anomalies which go beyond the scope of this paper. The justification of 
this approach is based on the fact that our considerations serve as a bridge between the manifestly 
noncovariant Dirac quantization and its Lagrangian path integral counterpart that can be cast into 
manifestly covariant form. It is the latter formalism that must be covariantly regulated to give a 
physically reasonable framework for rcnormalization and quantum anomalies. 
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the Pauli-Van Vleck-Morette formula |15[ - calculating the determinant of the matrix 
of second order derivatives of the principal Hamilton function with respect to its two 
arguments. However, in contrast with non-gauge theories this factor cannot be directly 
constructed in terms of this determinant, because in view of eqs.( |2.2j ) and ( |2.5| ) the 
matrix 

d 2 S(q,q>) 



Hk' 



dq % dq k ' 



has left and right zero-eigenvalue eigenvectors [17, 18 



0. 

= 0, 



dp k 



(2.6) 

and, therefore, is degenerate 

(2.7) 
(2.8) 



p = —dS/dq' 

The construction of the one-loop preexponential factor in terms of this degenerate 
matrix is equivalent to the Faddeev-Popov gauge-fixing procedure for gauge theories. 
It consists in adding the "gauge-breaking" term bilinear in "gauge conditions" Xf and 
Xy - two sets of arbitrary covectors at the configuration space points q and q' 



Sik' + X f C^Xy. 



(2.9) 



This allows one to replace the degenerate matrix S^i by the new invertible matrix D^t, 
provided that the gauge-fixing matrix C^ u is also invertible and these gauge covectors 



produce invertible "Faddeev-Popov operators" |2(J 



J = det J» ^ 0, 
J' = det J'^ ^ 0. 



In terms of these objects the solution of the continuity equations 
the following expression |T7|, [L^. |J 

n l/2 



det D 



J J' det C. 



flU _ 



(2.10) 
(2.11) 

is given by 



(2.12) 



which can be regarded as an analogue of the one-loop expression for the effective action 
of gauge field theory - the contribution of gauge fields det D ik i partly compensated by 
the contribution of ghosts J and J'. This compensation makes the prefactor indepen- 
dent of the introduced arbitrary elements of gauge- fixing procedure (Xf , X%,, C^ v ) - the 
analogue of on-shell gauge independence in gauge field theory. The mechanism of this 
gauge independence is based on the "Ward identities" for the gauge field "propagator" 

(U) 



Xy D 



-1 k'i 



(2.13) 



easily obtained by contracting (|2.9| ) with and using (|2.7|) . The use of these identities 
shows that arbitrary variations of the quantities (Xf , X k ,, C^ u ) in the one-loop prefactor 
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vanish due to the cancellation of terms coming from the gauge field det D^i and ghost 
J J' determinants. 

The nature of covectors (X^X^,) as matrices of gauge conditions does not only 
follow from the fact that they remove the degeneracy of the matrix So-/ caused by gauge 
invariance of the theory. As shown in [|T7|, |H| the quantum Hamiltonian reduction of the 
kernel K(q, q') with the one- loop factor ( |2.12| ) leads to the unitary evolution operator 
in the physical sector defined by the unitary gauge conditions 

= 0, (2.14) 

such that 

dX^ 

V - W - (2-15) 

The unitary gauge conditions are imposed only on phase space variables of the theory 
(q, p) (in this case only on coordinates In contrast with relativistic gauges involving 
Lagrange multiplyers, they manifestly incorporate unitarity and do not give rise to 
propagating ghosts. The price one usually pays for manifest unitarity is the absence 
of manifest covariance, that can be restored by going over to relativistic gauges. Thus, 



the solution ( 2.12 ) is obtained by directly solving the quantum Dirac constraints within 
the framework of unitary gauge conditions. In what follows we show that the same 
solution can be obtained by a direct calculation of the path integral in the relativistic 
gauge, and this derivation as a byproduct will establish the relation between unitary 
and relativistic gauge conditions. 



3. Lagrangian versus canonical formalisms 

As a first step towards the covariant path integral in the relativistic gauge let us 
consider the Lagrangian formalism of the theory with the canonical action For 
this we introduce the collective notation for the full set of Lagrangian configuration 
space variables - canonical coordinates and Lagrange multiplyers^ 

g a = (q*(t),N»(t)). (3.1) 

In what follows we shall need also (space)time condensed DeWitt notations in which 
the index a includes not only the spin labels and spatial coordinates x but also the 
time variable t, and the contraction of these indices implies the time integration (as 

2 Explicit time dependence of unitary gauge conditions is necessary in the theories with 
parametrized time in order to have nontrivial time evolution with a nonvanishing physical Hamil- 
tonian §, [[J. 

3 In the context of Einstein gravity theory this collection of ten fields g a ~ g a p(x,t) comprises 
the whole set of spacetime metric coefficients taken in a special parametrization adjusted to (3+1)- 
splitting: q l = g ab (x.,t), N a ~ g 0a (x.,t). 
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mentioned in Introduction we shall call them covariant) . In these notations the action 
has the form 

S[g]= dtL(q,q,N) (3.2) 

with the Lagrangian which does not involve time derivatives of the Lagrange multiplyers 
N = N^lt). This Lagrangian is related to the integrand of the canonical action ( |1 . 1| ) 

L{ q A^N)={ P ^-N^P))\ MN) (3-3) 

e substitutioi 
of the velocities q 



p=p°(q,q,N) 

by the substitution of the expression for the canonical momentum p®(q, q, N) in terms 



0/ • ,n dL(q,q } N) 
Pi (?, Q, N) = ^- (3.4) 

which is a solution of the canonical equation of motion 

dpi 

The classical action is invariant under gauge transformations with local (arbitrary 
time and space dependent) parameters / M = f^{t) vanishing on spacetime boundary 
f ll ('t±) — 0. For the canonical action ( II . 1| ) these transformations are canonical and, 
therefore, ultralocal in time for phase space variables, but involve the time derivative 
of the gauge parameter for Lagrange multiplyers 

5q l = {q\ W, S Pl = { Pl , T,}r, (3.6) 
SN" = f li - U£ u N a f u . (3.7) 

Thus the first class constraints T M (g,p) serve as generators of infinitesimal gauge trans- 
formations on the phase space of (q,p). Note, in particular, that the transformations of 
phase space coordinates q when restricted to the Lagrangian surface of the Hamilton- 
Jacobi function, p = dS/dq, read 

Sq* = V;/^ (3.8) 

in terms of the vector field ( |2.4| ) . 

These gauge invariance transformations in the Lagrangian formalism take the form 
of the infinitesimal transformations of configuration space variables 

6g a = R;r, (3.9) 

K 6 4^ = 0. (3.10) 
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with the Lagrangian generators Eft = (R^, R^) . Here we use covariant condensed 
notations in which the index a includes time, and its contraction implies the time inte- 
gration. Thus, i?" forms a delta-function type kernel with two entries a — > (a, t) , [i — > 
GM') 



pa 



R a Jd/dt) S(t - 1') 



(3.11) 



where R a (d/dt) denotes the differential (or ultralocal multiplication) operator acting 
on the first argument of the delta function. In view of the relation (|3.3| ) between the 
canonical and Lagrangian formalisms different components of this kernel follow from 
the invariance transformations ( |3.6|) - (|3.71) HQ 



r: 



dT 

dpi 




(3.12) 
(3.13) 



The distinguished role of the Lagrange multiplyers here manifests itself in the fact that 
only the a = a component of ( |3.11| ) forms the first order differential operator while the 
other components are ultralocal in time. In particular, the transformation (|3.12| ) is a 



Lagrangian form of (|3.8| ). 

In what follows we shall have to use condensed notations of both canonical and 
covariant nature. For brevity we shall not introduce special labels to distinguish be- 
tween them. As a rule, when the time argument is explicitly written we shall imply 
that the corresponding condensed index or indices are canonical, i.e. they contain only 
spin labels and spatial coordinates and their contraction does not involve implicit time 
integrals. For example, the left-hand side of gauge identities ( |3.10| ) can be written 
down in the form 



R 



a SS[g] 
5g a 



K{d/dt) 



S_S[g] 
6g*(t) ' 



(3.14) 



where the time integral implicit in the contraction of the covariant condensed index a 
removed the delta function contained in R^ and, thus, the result boiled down to the 
action of the differential operator R£(d/dt) on SS [g]/5g a (t). This operator obviously 



differs from that of eq. (|3.11| ) by the functional transpositon - integration by parts, 
because in contrast with ( |3.11[ ) it acts on test function with respect to upper condensed 
index a. This fact is indicated by the order of operator indices reversed relative to eq. 

(ED- 

4 The equivalence of the full set of transformations (137 



)-( |3.7| ) (including transformations of mo- 
menta) to (3.9)-(3.13) holds only up to terms proportional to equations of motion ||. Only up to 
such terms holds the equality Sp \ p °(q,q,N) = $P For us, however, it is only important to 

know that Lagrangian gauge transformations can be obtained from the subset of transformations of 
coordinates q l and Lagrange multiplyers N 11 in canonical theory. 
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Another important distinction between these two types of condensed notations con- 
cerns functional derivatives. We shall always reserve functional variational notations 
5/5g a = S/6g a (t) for variational derivatives with respect to functions of time, while 
the variational derivatives with respect to functions of spatial coordinates will be 
denoted by partial derivatives. For example, in the gravitational context we have 
S/5g a = 5/5gap(x,t) vs d/dq i = 5/5g ab (x). 



4. The path integral in relativistic gauge 



The path integral solution of quantum Dirac constraints was constructed in [11 



This is a functional integral over the gauge fields g a and grassman ghost fields C = C M 
and C = C v with a conventional Faddeev-Popov gauge fixing procedure 

K(q + ,q-)= J DgfjL[g]DCDCexp^(s[g]-^c^+C li Q^C v Y (4.1) 

The total action in the exponential here contains the gauge fixed classical and ghost 
actions. The gauge fixed classical action 

Sg[g] = S[g]-~x ti c^ X u , (4.2) 
includes the gauge breaking term 

I x'W = I f + dt x^g, g) c^x u (g, g) (4.3) 



2 ' v p " 2 
quadratic in the relativistic gauge conditions 

X" = X"(0,<7), (4-4) 
c£ = deta^ ± 0. (4.5) 

An important distinction of these gauge conditions is that they involve the velocities 
of all Lagrange multiplyers which results in propagating nature of all gauge and ghost 
fields. The corresponding Faddeev-Popov operator Q£, the kernel of the ghost action 
bilinear in ghost fields C and C 

QZ = %;K, (4.6) 

turns out to be a second order differential operator (cf. eqs.( |3.13| ) and (|4.5|) ) 

Q" u {d/dt)6{t - t') = f-< d 2 /dt 2 + ...) 6(t - t'). (4.7) 



An important ingredient of the path integral ( |4.1|) is the local integration measure 
l~t[g] which is built of the Hessian matrix of the classical Lagrangian and invertible 
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matrix c^ v fixing the gauge in ( |4.3|) 

H[g\ =l[[detG ik {t) detc^(t)] 1/2 = ( Det G ik Det c>) 1/2 , (4.8) 
t 



Note that the symbol det in ( |4.5| ) and in the middle part of eq.( [4.8[ ) denotes the 
determinants of matrices acting in the space of canonical condensed indices. The 
products over time points of local factors det Gi k (t) and det c^ u (t) can be regarded as 
determinants of higher functional dimensionality if we redefine these matrices as time- 
ultralocal operators acting in the space of condensed covariant indices Gik = Gik8(t—t') 
and c^y = Cp V 8(t — t'). We shall denote such functional determinants for both ultralocal 
and differential operators in time by Det. Thus the right-hand side of eq. ( |4.7|) is 
formally written down in terms of such (purely divergent) functional determinant of 
the ultralocal operator. For the Hessian matrix it looks like 

Det G ik 5(t - t') = exp | 5(0) J* + dt In det G ik J . (4.10) 

Similar definition holds for the ultralocal gauge-fixing matrix^. 

The final igredient which accomplishes the definition of the path integral fl4.1|) is a 
specification of boundary conditions on integration variables. Integration in ( |4.1|) runs 
over field histories with fixed values of canonical coordinates and ghost fields at t± 

q'(t±) = qi, C"(t±) = 0, C v (t ± ) = 0, (4.11) 

the coordinates q± being the arguments of the two-point kernel K(q + ,q_). On the 
contrary, the boundary values of Lagrange multiplyers are integrated over in the infinite 
range 

- oo < N»(t±) < +oo. (4.12) 

By using these boundary conditions, invariant with respect to BRST transforma- 
tions of the total action in ( |4.1| ), one can show the gauge independence of the path 
integral and also give a formal proof that this integral solves quantum Dirac con- 
straints in the coordinate representation of canonical commutation relations for q = q+ 
and p + = hd/idq + . This proof is based on an obvious consequence of the integration 
range for A^ M (t + ) 

I DsDCDC, mkr)^>°' (413) 



5 We don't consider here the case when is a differential or nonlocal operator. The only 
difference in this case is that it is no longer a part of the local measure. It contributes to the path 
integral a nontrivial functional determinant usually represented as a gaussian integral over auxiliary 
Nielsen-Kallosh ghosts. 
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where ellipses denote the full integrand of the path integral (|4.1|) . The functional 
differentiation here boils down in the main to the deexponentiation of the constraint 

^ = - W 0, P (* + ))U W) (4.14) 

(the differentiation of gauge-breaking and ghost terms gives terms cancelling one an- 
other in virtue of Ward identities []TT|). This deexponentiated constraint can be ex- 
tracted from under the integral sign in (|4.13|) in the form of the differential operator 



acting on q + , T M (g + , hd/idq + ), so that the equation (j4.13|) takes the form of the quan- 



tum Dirac constraint |14[]. This nonperturbative derivation is, however, purely formal 



and in an uncontrolable way depends on the skeletonization of the path integral |I4 
In contrast with these formal considerations in the following sections we prove that 
this path integral solves the quantum constraints at least in the one-loop order of 
semiclassical expansion. 



The Feynman diagrammatic technique for the integral ( |4.1| ) was also built in [|TI 
Its main emphasis concerns, certainly, the boundary conditions at t±, because in other 
respects the ^-expansion produces a standard set of Feynman graphs. The first step 
in this expansion consists in finding the stationary point of the path integral subject 
to boundary conditions A4.11|) - ( |4.12j ). As shown in [|TTf| , for gauge fixed action (|4[ 



this is a unique solution g — g(t \ q + , g_) of the following boundary value problem for 
classical equations of motion in the chosen relativistic gaugef] 

T-^ri = 0, 4.15 

x ^g,g) = 0, t.<t<t+, (4.16) 
q(t±) = q ± . (4.17) 

Performing the gaussian integration over quantum fields in the vicinity of this sta- 
tionary point one arrives at the semiclassical answer (|2.1[) for our path integral with 
the tree-level action 

S{q,q') = S[g(t\q + ,q-)) (4.18) 

and formal expression for the one-loop preexponential factor in terms of the functional 
determinants of the ghost and gauge field operators and the local measure 

Det Q» 



V ao ) g=g(t\q + ,q-) 

6 The variation of the gauge fixed action gives for t- < t < £+ classical equations amended by the 
gauge fixing term and two conditions at the boundaries t = t± (requirement of vanishing coefficients 



of SN (1 (t±)). In view of (4.5) these conditions reduce to gauge conditions x 1 * — enforced at t = t±. 



On the other hand, in view of the identity ( 3.1 0| ) the classical equations with a gauge fixing term 



result in homogeneous equations for these gauge conditions Q^c^ a x a = 0, t_ < t < t + . The second 
order differential operator Q£ is assumed to be invertible under the Dirichlet boundary conditions, so 
that gauge conditions are enforced for all t. 
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The gauge field operator is given by the second order variational derivatives of the 
gauge- fixed action (|4.2|) . In the covariant condensed notations it equals 



\b = S ah - X >1 aC,uXl (4-20) 
S 2 S[g[ 
5g a 5g b 



Sa b = T377^ (4-21) 



- ^, (4.22) 

where the functional matrix \a °f linearized gauge conditions is a first-order differential 
operator with the delta-function type kernel 

X£ = xr a (d/dt)8(t-lf). (4.23) 

A formal definition of functional determinants in Q4.19Q is incomplete unless one 
specifies a functional space on which they are calculated. This is equivalent to speci- 
fying the boundary conditions for Green's functions of the gauge and ghost operators 
Qba _ p-iba an( j Q-iu w j 1 j c j 1 participate in the variational equations for their deter- 
minants [21] 



5 In Det F ab = G ba 5F ab , 5 In Det = Q' 1 V 8Q$ (4.24) 

and multiloop graphs of semiclassical expansion. 

Boundary conditions for the gauge propagator were derived in They actually 



follow from the linearized version of the boundary value problem (|4.15|) - (|4.17|) . The 
propagator G ab (t,t') is (n + m) x (n + m) matrix-valued function (cf. eq. (|1.8| )). Thus, 
it requires (n + m) boundary conditions at t±, n of them obviously being the Dirichlet 
conditions for its a = z-components (i = 1, ...n), while the rest of them coinciding with 
m linearized gauge conditions at t — t± 

F ca {d/dt)G ab {t, t') = 5 b c 5{t - t'), (4.25) 
G ib {t±,t') = 0, (4.26) 
X fl a (d/dt ± )G ab (t ± ,t')=0. (4.27) 



The latter belong to Robin type because for relativistic gauges the operator ( |4.23|) 



contains derivatives transversal to the boundary. These boundary conditions have the 
properties of the BRST-invariance, guarantee gauge independence of the path integral 
( |4.1| ) [jnj| and selfadjointness of the gauge field operator recently discussed in [ p2| . 



Finally, the propagator of ghost fields in ( J4.24 ) satisfies in view of ( [4.11] ) the Dirichlet 
boundary value problem 

Qf{d/dt)QZ 1 0(t,t) = tfl, Q- a ip (t ± ,t')=0. (4.28) 
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5. Local measure, Ward identities and gauge inde- 
pendence 

The prescription of boundary conditions for Green's functions does not, however, 
uniquely determine the variation of the determinants ( [4.24 ) . Problem is that kernels 



of propagators are not smooth functions of their arguments, and their irregularity 
enhances when they are acted upon by two derivatives contained in SF a b(d/dt) and 
5QJ 1 {d/dt). Therefore, one has to prescribe the way how these derivatives act on both 
arguments of Green's functions and how to take the coincidence limit of the resulting 
singular kernels in the functional traces of ( 4.24|) . It is well known that strongest 



divergences arising from these singular kernels are compensated by the local measure 



fj, [g] p3[. As far as it concerns the rest finite (or less divergent in field theories) part, 
its construction should maintain the gauge invariance properties encoded in the path 
integral at nonperturbative level. This might serve as a hint that can (at least partially) 
fix the ambiguities of the above type. For this reason, in this section we consider the 
role of local measure and Ward identities in the construction of the gauge-independent 



prefactor ( [4 . 1 9[ ) . 



It is well known that strongest power divergences of functional determinants are re- 



lated to terms with second order derivatives in the ghost ( [4.71) an d gauge field operators 



n in 



F ab (d/dt) = -^-a ab ^- + (5.1) 
at dt 



««* = t^t. (5.2) 



dg a dg l 

The Hessian matrix of the gauge-fixed Lagrangian is nondegenerate due to the rela- 
tivists nature of gauge conditions ( |4.4| )-(4~5). The gauge breaking term contributes 



the part quadratic in velocities of Lagrange multiplyers that were initially absent 
in the original Lagrangian L(q, q, N) 

L gf(g, g) = Hq, q, N ) - ^(g, g) c^x u (g, g)- (5.3) 

Therefore, the Hessian matrix has the following components 

W (5 ' 4) 

(5.5) 

(5.6) 

where is the Hessian matrix ( |4.5| ) of the ghost Lagrangian. In virtue of the relation 

det a a b = det Gik det c a p ( det ) 2 (5.7) 





— ' 






_ d X a 


c a /3 a J, 








= —a a 


c a/ 3 al, 
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the expression for the local measure Q4.S| ) can be rewritten as the following ratio of the 
determinants of the gauge field and ghost Hessian matrices 

Deta ab Y / 2 
Det at 

Therefore the one-loop prefactor (|4.19|) takes the form 



m^] = v • (5.8) 



_ / DetF afc \- 1/2 DetQ£ 



(5.9) 

9=9(t\q+,q-) 



This form is especially adjusted to the manifest cancellation of strongest divergent 
parts of determinants: as shown in []TB the <5(0)-type divergences for a second order 



differential operator are cancelled by the functional determinant of the corresponding 
ultralocal matrix coefficient of its second order derivatives. 

The representation ( |5.8|) for the local measure, which is similar to the ratio of gauge 
field and ghost functional determinants in ( [4.19| ), presents a simplest demonstration of 
gauge independence. The numerator and denominator separately depend on the choice 
of gauge conditions x^id^d) (the corresponding matrix = —dx^/dN"), but in the 
ratio this dependence cancels out and the total local measure ( |4.8| ) turns out to be 
gauge independent^ . The mechanism of gauge independence for the rest part of the 
one-loop prefactor is more complicated and is based on Ward identities for gauge and 
ghost propagators. 



These identities follow from the gauge invariance of the classical action (|3.10| ). 



The functional differentiation of fl3.10|) shows that on shell, that is on the background 



R;S ab = S a -f = (5.10) 



satisfying classical equations of motion, the functional matrix S ab is degenerate because 
it has zero-eigenvalue eigenvectors - the gauge generators 

SB?, 
Sg b 

(cf. the analogous property ( [2.7|) in the canonical context). As a consequence the gauge 
operator F ab satisfies the relation 

R»F ab = -Qlc af}X l (5.11) 

which can be functionally contracted with matrices of the gauge G bc and ghost Qa lfl 
propagators introduced above. Integration by parts of the derivatives in the operators 
= R a (d/dt) and = Q^(d/dt) does not produce additional surface terms in view 
of the boundary conditions for the propagators. Therefore, after using the equations 
C|4 . 2 5|) and (|4.28| ) for gauge and ghost propagators one arrives at the needed Ward 
identity 

c a( s X P b (d/dt)G bc (t,t') = -Q-^it^R^d/dt 1 ). (5.12) 



7 Gauge independent expression for local measure (denned by the Hessian of the classical gauge 
action (|4.9|)) was first obtained in the context of Einstein gravity theory m and then derived for 



generic gauge theory from the canonical path integral || . 
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Although this identity can be obtained from Q5.11|) by formal inversion of F ab and as 



finite-dimensional matrices, we emphasize the necessity of a cautious use of condensed 
notations with regard to possible surface terms following from integration by parts. The 
result is presented in the form clearly indicating the differential structure of operators 
acting on the arguments of Green's function. In what follows we shall often label the 
differential operators by arrows to show the direction (right or left) in which they act. 

The obtained Ward identities lead to gauge independence of the one-loop prefactor 
(|4.19|) provided that we consistently fix the functional composition laws for G ba 5F ab 
and Q~^ ll/ 5Q% in (|4.24j ). Let us assume that the time derivatives of varied differential 
operators 5F ab (d/dt) and 5Q£(d/dt) are understood as acting in two different ways 

5\nBetF ab = G ba 5F ab , (5.13) 

SlnDetQ^Q-^SQ/. (5.14) 

In contrast with the ghost operator, for which both of its derivatives are acting on 
one argument of the Green's function, eq.( |5.13j ) here implies a symmetric action of 
5F ab (d/dt) on both arguments of G ba (t,t') in the sense that 

L$ = \v a {t)F ab (d/dt)<p b {t) (5.15) 

represents the quadratic part of the gauge-fixed Lagrangian in perturbations of field 
variables (p a . It contains the squares of their velocities ip a (t) rather than their second 
derivatives (and SF ab (d/dt) obviously implies the variation of this operator with respect 
to its background field dependence). 

With these conventions the gauge variation of the gauge field and ghost determi- 
nants 

5 X In Det F ab = 5 x F ab G ba = -2 c Q/3 X I G ba 5x a a 



-2 J dt 



c aP xl{d/dt) G ba (t,t') 5X a a (d/dt')\ t=ti , (5.16) 

6 x \nDetQ» = Q^b x Q% = Q~ a ip R b ^X a b 

r Q- ip (t,t')R b p (d/dt')5xt(d/dt') 



dt 



(5.17) 



t=t' 



cancel out in virtue of the Ward identities (|5.12|) in the one-loop prefactor fl4.19|) which 



proves its gauge independence. Its independence from the choice of the gauge fixing 
matrix c^ u is also based on these identities, though in this case it is cancelled by the 
Nielsen-Kallosh factor in the local measure. 

For the one-loop effective action these Ward identities were in much detail analyzed 
24| in the framework of formal condensed notations. Here the main emphasis in the 



m 



mechanism of Ward identities is focused on boundary conditions and accurate definition 
of the functional determinants in ( [4.24|) . Another choice of the functional composition 



law in these variational equations leads in general to extra surface terms breaking the 
gauge independence of the one-loop prefactor. 
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6. Gauge field functional determinant 



Here we calculate the contribution of the gauge field determinant to the one-loop 
prefactor fl5.9Q . For this purpose we introduce matrix notations for operators acting in 
the vector space of canonical indices. In these notations the gauge field operator 



F{d/dt) = F ab (d/dt) = 
acts in the space of columns 
ip{t) 



F ik (d/dt) F iv {d/dt) 
F^d/dt) F^d/dt) 



Generically it has the form of the second order differential operator 

+ c, 



^ ^ ^ dt dt dt dt 



(6.1) 



(6.2) 



(6.3) 



where the coefficients a = a a b(t), b = b a b{t) and c = c a b(t) are the matrices acting in 
the vector space of <p (£), and the superscript T denotes their functional transposition 
(b T ) ab = bba- These coefficients can be easily expressed as mixed second-order deriva- 
tives of the Lagrangian in the gauge-fixed action ( |4.2j) with respect to g a and g a . In 
particular, the matrix of second order derivatives a a b is given by the Hessian matrix 



In accordance with the notation ( |5.15|) the "integrated by parts" version of the 
operator (|6.1|), F (d/dt), determines the quadratic part of the gauge fixed Lagrangian 
with two time derivatives of F acting symmetrically on two functions <p(t) and its 
transpose <p T {t). For arbitrary two test functions ipi and ip 2 the operator of the form 
( |6.3f ), in our notations, gives rise to the bilinear form 



ip l Fip 2 = 0! a 02 + 0i b<f2 + <fi b 02 + ¥>i cip 2 
and implies the following integration by parts 

<plF<p2 = ipl (F<p 2 ) + (W0 2 ) 
W = W(d/dt) = aj t + b. 
Here W is the Wronskian operator which enters the Wronskian relation for F 

ifl (F<p 2 ) - (Fipifipz = -— [ip[ (W<p 2 ) - (W Vl ) T V2 



(6.4) 

(6.5) 
(6.6) 

(6.7) 



and also participates in the variational equation for the canonical momentum dL^/dg 
valid for arbitrary field variations 5g (t) 

.<9L gf 



6- 



dg 



W(d/dt) 5g (t). 



(6.8) 
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Since the velocities of Lagrange multiplyers enter L g f(g,g) only through the gauge 
breaking term, the /i-component of the Wronskian operator is given by 

d X ° 



W^d/dt) 



c a/ 3 X b (d/dt), 



(6.9) 



where, in particular, the differential operator of the linearized gauge conditions Q4.23Q 
coincides with that of the boundary value problem ( |4.27p . 

Now we introduce the notation for the matrix valued Green's function of this bound- 
ary value problem 



G ab (t,t') 



G ik {t,t') G iu {t,t') 



= G(t,t'), 



(6.10) 



and in accordance with equations fl5.13| ) and (|6T4| ) write down the variation of the gauge 
field determinant 



cHnDetF = Tr 5FG 



dttr 



dt ti 



d 



5FG(t,t) 



d 



(5a -V + 5b -f - + 5b T 4- + 5c) G (t, t) 

y dtdt' dt' dt ' v ' ; 



(6.11) 



t'=t 



Here tr denotes the matrix trace operation with respect to condensed canonical indices 
of 5a = 5a ab (t), 5b = 5b ab (t), 5b T = 8b ba (t), 5c = 5c ab (t) and G (t,f) = G ab (*,*')• 

Our goal now will be to construct a special representation of the Green's function 
G(t,t') and integrate the variational equation ( |6.11| ) by the method of JT6[. For this 
purpose, similarly to [|nj, we introduce two sets of basis functions u_ and u + of the 
operator F 



Fu ± = 0, u± = u± a A (t) 



(6.12) 



satisfying the full set of boundary conditions ( |4.26|) - ([4.27|) respectively at t- and t + . 
In view of ( |6.9[ ) and invertibility of the matrix dx a /dN' 1 the Robin type boundary 
conditions ( |4.27| ) can be rewritten in terms of the /x-components of the Wronskian 
operator, so that the full set of these boundary conditions takes the form 



U + \(t + ) =0, W^aU+ a A (t + ) =0, 
uJ A {tJ) =0, W^aU- A (t_) =0. 



(6.13) 
(6.14) 



In condensed notations we regard these basis functions, enumerated by the condensed 
index A of arbitrary nature, as forming the square matrices with the first (contravari- 
ant) index i and the second (covariant) index A. 

The Wronskian relation ( |6.7| ) can be used to form the t-independent matrix of 
the Wronskian inner products A i2 = (pj (Wy3 2 ) — (Wiy9i) T ^2 of basis functions <pi t 2 = 
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u±(t). In view of the boundary conditions of the above type this matrix has only two 
nonvanishing blocks given by the two mutually transposed matrices 



A+- 



u 



T 

+ 

T <Wu A 



U 



:wu_) T u, 



(A-+)ab, 



This fact can also be represented in the form of the following matrix relation 



u 



u 



T 



[Wu- 
(Wu A 





Wu + WuJ 




"A_ + 




u + u 




A + __ 



(6.15) 
(6.16) 
(6.17) 



(6.18) 



Under the assumption of absence of zero modes of F subject to boundary conditions 
fl4.26p -( f07D (absence of linear dependence of u + and u_) this relation implies the 
invertibility of A_ + . It also allows one to to establish the following important relations 
for equal-time bilinear combinations of basis functions |l6j 



u+(t) (A_ + )-^ T (t) + u_(f) (A+^uKt) = 0, 
u + (t) (A_ + )~ V(t) + «_(t) (A + _y l u T Jt) 



a 



(6.19) 
(6.20) 



Here 1 = 5% denotes the unity matrix in the space of canonical indices a. To clarify 
their use in these equations and in what follows, we note that in the transposed matrix 
= u± a A the covariant index A is considered to be the first one (in contrast with 



-1 TAB 



gives rise m 



u±), so that the matrix composition law with (A + _) = [ (A_ 
eqs.( 6.20 ) to the matrices with two indices a and b. 

By the method of ref . fl6|] (applied there in the case of Dirichlet boundary conditions) 
one can show that the Green's function of the mixed Dirichlet-Robin boundary value 
problem ( |4.25| )-( [i.27| ) has the following representation 



G(t,t') = -6 (t-t')u+(t) (A_ + ) u_(t' 



+e (t'-t)u4t) (A+^uKt') 



(6.21) 



where 9{x) is the step function: 9(x) = 1 for x > and 9(x) =0 for x < 0. This ex- 
pression is the analogue of positive-negative frequency decomposition for the Feynman 
propagator in scattering theory. 

Substituting ( |6.21| ) to ( |6.11[ ) and repeating the calculations of |TB[ one can see that 
in view of relations ( |6.19|) -( |O0| ) the 5(0)-type terms with derivatives of step functions 
reduce to the variation of the logarithm of the local measure for F. Therefore, the 
variational equation Q6.ll ) acquires the form taking into account the cancellation of 
<5(0)-type divergences 



5 In 



DetF 
Det a 



-tr (A. 



dt u T 5F «j 



(6.22) 
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Here we reserve the same notation tr for the trace operation in the space of indices A 
enumerating the basis functions of F, so that the matrix multiplication here should 
read [ (A_+)- x } AB u T B 5F u +A . 

The further transformation of the time integral in ( ft-22 ) is based on integration by 
parts which uses the identities (|6.5|) , (|6.7|) and their corollaries obtained by replacing the 
operators F and W with their variations 5F and 5W. This allows one in a systematic 
way to reduce the integrand to a total derivative modulo the terms vanishing in virtue 
of equations (|6.12| ) for basis functions or their varied versions 5Fu± = —F5u±. The 
result boils down to the contribution of surface terms at t = t±. In view of boundary 
conditions for u± they takes the following form 



dt u T 5F u A 



ul5(Wu + )i - 
+ \(Wu. 



iSu\ - ul5(Wu 



where for brevity we introduced notations 

(Wu±)i = WiaU±, {WU±) 11 = W„ a U ± . 



(6.23) 



(6.24) 



It is useful to rewrite this result in the matrix form 



dt u 5F u A 



u 



5(Wu. 



6(Wu + ), 



(6.25) 



For brevity here the indices A,B, ... enumerating the basis functions are omitted: they 
can be regarded as included into subscripts ± - the rule which we shall imply in what 
follows. Thus all the matrices in this relation are square (but certainly not symmetric). 
One index of these matrices, with respect to which the contraction takes place, is a 
combination of i and /i (being in superscript and subscript positions respectively or 
vicy versa), while another index is A encoded as we agreed in ±. 

The expression for the time- independent Wronskian matrix A_ + can also be rewrit- 
ten in terms of similar matrices. When calculated at t + and t_ it looks respectively 
as 



u 



[Wu-X ut 



{Wu+)i 



u 

Wu 



(6.26) 
(6.27) 



j t- 



The nondegeneracy of A_ + guarantees the invertibility of matrix factors in these two 
expressions. 
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Combining the obtained expressions (|6.25|) , (|6.26|) , (|6.27|) with equation (|6.22|) one 
can convert its right-hand side into the total variation 



5 In 



DetF 
Det a 



-5 tr In 



'{wu+y 


+ 5 tr In 













S In det D (6.28) 



of the logarithm of the determinant of the matrix D = with indices a 
b = (k', v) 



D 



Wu 



[WuJ)y Ul 



(i, //) and 
(6.29) 



D 



The block form of this matrix 
Sik' X" 

Sw = -(Wu+UU) All (^_),(t) 
X£ = <(t+)A:i(Wu_)*(t_), 

c«' = <(t + )A:i 

allows one to rewrite its determinant in the form 



det D 



det (Sw+XtC^XZ) detC^ 



C„ 



(6.30) 
(6.31) 
(6.32) 
(6.33) 

(6.34) 



The notation Siy chosen for the block ( |6.30|) of the matrix D is not accidental. This 
expression really yields the matrix of second-order derivatives of the Hamilton- Jacobi 
function ( |2.6j ) with respect to its arguments q l and q k ' . Indeed, this matrix coincides 
with the derivative of the canonical momentum — = dS/dq k ' with respect to q l . 
The momentum is taken at the initial moment t_ on the classical extremal joining the 
points q and q'. In view of the relation (|6.8j ) 

dg v (i_ | q+,q-) 



Sik' 



■W k 'b'(d/dt. 



dq\ 



(6.35) 



where the derivative dg b (t\q + , q_)/ dq\ is given in terms of the Green's function G ab (t, t') 
(see refJ^Bl) 



dg b (t\q + ,q_ 
dq\ 



-W ia (d/dt + )G ab (t +} t) 



(6.36) 



whence one gets the equation ( |6.30| ) after using the basis functions representation for 
G ab (t + ,t). 

Thus, up to field independent normalization constant, the one- loop contribution of 
gauge fields equals 



DetF 
Det a 



Const 



det ( S lk , + XtC^Xl, 



det C, 



(6.37) 
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This is a part of expression (|2.12| ) for the one-loop prefactor provided we identify arbi- 
trary elements (Xf , X%,, C^ u ) of its canonical gauge fixing procedure with the quantities 
(|6.31|) , ( |6.32|) and (|6.33|) above. Let us now show that the contribution of ghost fields 
gives the remaining part of ( |2.12| ). 



7. Functional determinant for ghost fields 



The ghost operator (|4.6| ) is not even formally symmetric because its right and left 
action is defined correspondingly on contravariant and covariant vector fields in the 
space of gauge indices 

Q^(d/dt)F(t) = X\{d/dt) R a v {d/dt) /"(*), (7.1) 
Q/id/dt) f,(t) = R u a {d/dt) x£(d/dt) U(t), (7.2) 

where the form of the first-order differential operators R^d/dt) and x%{d/dt) varies 
depending they are acting on test functions (<p, f) with respect to their condensed 
indices aor/j 

X \{d/dt) <f(t) = + ^(t), (7.3) 



dg a dt dg a 

d dx 11 d\ 
dt dg a dg a 



X/(d/dt) f,(t) = I + I Ut) (7.4) 



and 



R a jd/dt)nt) = (s a ± + ..) nt), (7.5) 



dt 



R u a (d/dt) <p a (t) = (-5lj t + Mt). (7.6) 

Pairs of relations ([7.1|)-( |7T2"D , (|7.3|)-( f7T4|) and (|7.5|)-fl7lf) obviously differ from one an- 
other by integration by parts in bilinear integral forms with Q^, Xa an d Rt as functional 
two-point kernels. 

Because of the absence of symmetry the Wronskian relation for 

d 



fiAQ"J2)-(fi,Q^)f2=~ dt 



(7.7) 



involves two different Wronskian operators W^ v and W^ u with the arrows indicating 
the direction in which the differential operators are acting. Their actions on test 
functions 

w\r = R\(d/dt)nt) (7.8) 

U ™\ = W/U = ( + ^ J U(f) =X/{d/dt) f^t) (7.9) 
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are essentially different from one another and do not differ by a simple functional 
transposition of one and the same operator^]. Note, in particular, that the operator 
(|7.9|) coincides with with the ^-component of the linearized gauge operator ( |7.4|) . 

With these notations one can again construct the Green's function of with 
Dirichlet boundary conditions in terms of doubled set of the right and left basis functions 
r£(*) and v±(t) 



QV±(*) = o, r£(t±) = 0, 
Q l ?v±{t) = Q, <(t±) = 0. 



(7.10) 
(7.11) 



The matrix of Wronskian inner products of these basis functions has only the following 
nonvanishing block components 



e+ = rt (w;v+)-(rtw;)v+, 
o+=r$ {w;%)-Kw;)v- 

participating in the analogue of matrix relations ( |6.18| ) 



■{Wr. 
-(Wr, 



(Wv+)^ (WV 





~6+ 




9 + ~ _ 



where for brevity we introduced the notations (Wr±) v = r± W 



(7.12) 
(7.13) 

(7.14) 



W v n r ± and 



{yVv^^ = W u v v ■ Similarly to (|6.19| ) these equations imply the equal-time bilinear 
relations for ghost basis functions 

{Wv + ) IM (9_ + )- 1 r1 + (Wv-) fi (9 + ~ 
«+ (6+y 1 (Wr-T + v- (Wr, 

vt (9+y 1 r u _ + v~ (err 1 r" = 0. 



v 



In terms of these basis functions the ghost Green's function equals 



(7.15) 
(7.16) 
(7.17) 



Q- lv (t,t') = -9(t-t')vZ(t) (9_ + r 1 rl(t') 



+e(t'-t)v-(t) (err'rUt'). 



(7.18) 



Substituting this expression to ( |5.14j ) and again repeating the calculations of |16|] one 
gets the <5(0)-type terms absorbed by the local measure and arrives at the expression 



5 In 



Pet Q£ 
Det a„ 



tr 



dtr"_(t) 5Q/v+(t) 



1 

+ 2 tr ' 



t- 



dtrl(t) SQ/v~{t), 



(7.19) 



8 We shall not introduce for these two Wronskian operators separate notations and will only 
distinguish them by the (upper or lower) position of their indices. The order of indices depends on the 
direction in which the operator is acting: when acting to the right its second index (irrespective of its 
covariant or contravariant nature) gets contracted with the index of the test function. For operators 
acting to the left the order of indices reverses, so that with these conventions W^ v =W J 1 '■ 
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where the 1/2 coefficients originate from the symmetric prescription for the theta func- 
tion| 9(0) = 1/2. It differs from the analogous expression ( 6.22j) by the presence of two 
terms built of two different sets of right and left basis functions - the consequence of 
asymmetry for the ghost operator as compared to the symmetric gauge field one. These 
terms are different here in contrast with the case of for which they coincide and, 
thus, lead to 1/2 coefficients adding up to unity (or, irrespective of the rule #(0) = 1/2, 
lead to using the identity 9(t — t') +9(t' — t) = 1). Their further transformation repeats 
the calculations of fl6j briefly explained in the previous section and leads to the result 



-tr 



tr i 



dtr*_(t) 5Q/v+(t) 



5 In det 



+ (t_)(^+)-V^(t + ) +< yindet<(t 



dtrl(t) 8Q v »v-{t) 
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(7.20) 



(7.21) 



Each of the first terms on the right hand sides of these two relations is equivalent to 



the equations (3.19)-(3.20) of the reference [16[] for a symmetric operator subject to 
Dirichlet boundary conditions. On the contrary, the terms with det a^(t±) are new 

and originate from asymmetric action of 5Q J 1 in the variational definition of the ghost 
determinant ( [5. 14| ) . Using 9_ + and 9 + ~ calculated from ( |7.12| )-( [7.13| ) respectively at t- 
and t + and taking into account that (Wr±) fl (t±) = a£ (t±)r u (t±) one finds the following 
expressions for matrices 



<(*-)(0- + rv(t + ) 



-ir-(U) 



(*-)]; 
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(*-)] 
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(7.22) 
(7.23) 



Here we imply that r±(t) = r±(t) and r±(t) = f±(t) are square matrices with the 
first index \i and the second index encoded in their subscripts ±. Substituting these 
expressions to (|7.20|) - (|7.21|) and then to (|7.19| ) one finds that the variations of det a^(t±) 
cancel out and the result can be functionally integrated 



DetQ£ 



const ( det J% det J? ) 



-1/2 



Det at 

JZ=[r-{t-){r-y\t + )Y v , 

^ = -[^ + (M(r + )- 1 (t-)]: 



(7.24) 

(7.25) 
(7.26) 



up to numerical normalization constant. The coincidence of notations for the last two 
matrices with those of Faddeev-Popov operators in unitary gauges ( |2.10|) -( p.ll| ) is again 
not accidental. The proof of these equalities is as follows. 

9 This prescription is apparently related to symmetric ordering of equal-time operator products 
which makes them Hermitian, but we shall not trace back to the roots of this rule in operator 
quantization. 
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As it follows from the previous section, to compare our one-loop prexponential 
factor with that of the solution ( 2.12 ) we have to identify the quantities ( |6.31 )-( |6.32 ) 
with matrices of canonical gauges. The Faddev-Popov matrix corresponding to the 
gauge matrix ( |6.31| ) is therefore 

J£ = X? Vt = V* (Wu + Ut + ) AI^(*_). (7.27) 

We show now that this expression equals ([7.25 ). For this purpose, note that the 
combination VJ, (Wu+)i(t + ) can be identified with the linearized constraint under the 
variation of canonical coordinates and momenta induced by variations of Lagrangian 
variables 5g a (t) = 

8q i (t + )=u\(t + ) = 0, (7.28) 
5 Pi (t+) = Wf a (d/dt + ) u a + (t) = W ia (d/dt+) u\(t). (7.29) 

Here Wf a (d/dt) is a Wronskian operator built similarly to ( p.8| ) but with respect to 
the original Lagrangian L(q, q, N) instead of the gauge-fixed one. The second equality 
in ( [7.291) follows from the Robin-boundary conditions on basis functions (|6.13|) . Thus, 
since is the momentum derivative of the constraint (|2.4| ), the combination of the 
above type coincides with the linearized constraint which in view of (|4.14j) equals the 
variation of the /i-component of the classical equations of motion 

Vt (Wu+UU) = STMt+lpiU)) = -S, a (d/dt + ) u a (t + ). (7.30) 

Then, in virtue of equations ( |6.12| ) for basis functions and relation (|7.9|) for the Wron- 
skian operator of Q^(d/dt) 

Vt(W«+) i (t + ) = - W u a (d/dt + )(c aP X^ a u a + )(t + ). (7.31) 

Substitution of this relation to (|7.27|) leads to the expression bilinear in basis functions 
u±(t) with one of them acted upon by the operator of linearized gauge conditions. Such 
an expression can be transformed by using the Ward identitiy ( |5.12|) . Substituting the 
basis function representations ( |6.21|) and ( |7.18| ) into it at t > t' one gets 

c a(3 xP(d/dt)u a + (t) a:X-(0 = -v+{t) {e+y'r^R^d/df). (7.32) 

Using this relation with t = t + and t' = t_ together with (|7.31|) in the right hand 
side of ( |7.27|) and using the Dirichlet boundary conditions for ghost basis functions one 
finally arrives at the equation ( |7.25| ) . Similar proof holds for ( |7.26|) . 

This accomplishes the derivation of the equation (|7.24|) . It represents a reduction 
algorithm that expresses the ghost functional determinant in terms of determinants of 
the canonical Faddeev-Popov matrices with qualitatively lower functional dimension- 
ality. Simultaneously this algorithm relates the gauge fixing procedure in relativistic 
gauges to the special unitary gauges of the form ( |2.14D with the matrix ( |2.15| ) given 
by ( 6.31|) . This reduction algorithm together with the the result for the gauge field 
determinant (|6.37j) finalizes the the proof that the path integral with the one-loop pref- 
actor ( |4.19| ) indeed represents the semiclassical solution (|2.1| ), (|2.12| ) of quantum Dirac 
constraints. 
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8. Conclusions 



The equality of expressions (|2.12|) and ( }4.19|) comprises the equality of two gauge 
independent objects - the function ( |2.12| ) independent of the unitary gauge conditions 
(Xf ,X%) and the functional ( }4.19| ) independent of the relativistic ones \a- The mean- 
ing of this seemingly vacuous equation is that it establishes the mechanism of identical 
transition from unitary to relativistic gauge conditions. Such a transition is very im- 
portant because it proves intrinsic unitarity of a manifestly covariant quantization in 
terms of Lagrangian path integral. This transition is generally a very nontrivial proce- 
dure, because in contrast with unitary gauges in relativistic gauges the theory possesses 
qualitatively different number of propagating degrees of freedom: gauge modes of g a , 
Lagrange multiplyers as well as ghost fields (C 1 ', C v ) are dynamical. It is a subtle 
mechanism of BRST invariance that guarantees the cancellation of contributions of all 
these modes in physical quantities. This mechanism makes the result equivalent to 
quantization in the physical sector, that is in unitary gauge. However, this transition 
is usually performed at the level of formal identical transformation in the path integral, 
which is reached only in a singular limit e — > of the so-called e-procedure of ref. ||. 

On the contrary, the equivalence of ( [2.12|) and fl4.19|) is actually achieved without 
any singular limiting operations - by independently calculating these quantities in 
unitary and relativistic gauges and observing their coincidence when these two sets of 
gauge conditions are related by equations ( |6.31| )-( |6733D . These relations are nonlocal 
in time - the matrices of unitary gauges and gauge fixing matrix C^ v express in terms 
of basis functions of the gauge field operator, nonlocally depending on its relativistic 
gauge fixing elements \a an d c^ u . 

In physical applications, loop expansion of the path integral as a means of solving 
quantum Dirac constraints has a very important advantage. This technique implies 
that solutions of canonical constraints (attached to time foliation) admit the spacetime 
covariant description in terms of Feynman diagrams [p5|| . This property is of crucial 
importance for correct regularization of inevitable ultraviolet divergences. This regu- 
larization should maintain the covariance in the form not splitted by time foliation. 
Spacetime covariance is not manifest due to the canonical origin of Dirac constraints, 
but it gets restored in the proposed calculation technique, because the arising one-loop 
functional determinants can be cast into spacetime covariant form by a suitable choice 
of relativistic gauge conditions. Important implication of this technique is the theory 
of loop effects in quantum cosmology including, in particular, effective equations for 
expectation values in the early inflationary universe [^6] . 

The aspects of gauge independence considered above are important in gauge theory 
applications in spacetimes with boundaries or nontrivial time foliations. There exists 
a long list of examples when the calculations of a formally gauge independent quantity 
- one- loop effective action - gives different results in different gauges ]37], |2"8fl . No 
exhaustive explanation for these discrepancies has thus far been given, and there is a 
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hope that a careful analysis of Ward identities with regard to boundary conditions can 
resolve this problem. 

The last but not the least problem that belongs to the scope of our result is the 
theory of quantum gravitational tunneling and the physics of wormholes |[29|| . There 
is a widespread opinion that the predictions of Euclidean quantum gravity modelling 
these phenomena have a questionable status due to the indefmiteness of the Euclidean 
gravitational action [ 30f1 . The conclusion drawn in J3l| about the mechanism of tran- 
sitions with changing spacetime topology is based on the existence of a negative mode 
on the wormhole instanton - a formal extrapolation of the mechanism which is di- 



rectly applicable only to non-gravitational systems [32]. As it was pointed out in [31 



(see also |HJ), this issue should be revised from the viewpoint of the Wheeler-DeWitt 
equation. The necessity of such a revision is obvious because this negative mode be- 
longs to the conformal sector which is not dynamically independent in the Lorentzian 
quantum gravity. Its contribution is, therefore, cancelled by ghost fields in relativistic 
gauges. One should expect that similar cancellation should take place also in Euclidean 
theory (despite the hyperbolic vs elliptic nature of field equations in Lorentzian and 
Euclidean theories, they demand, after all, the same total number of boundary con- 
ditions)^ Other difficult ies in Euclidean gravity theory related to this issue also find 
place in current literature: the lack of strong ellipticity of the Dirichlet-Robin boundary 
value problem ( 4.25 )-( 4.27 ) observed in [34[ seems to be explained by the indefiniteness 
of the Euclidean gravitational action. The technique proposed here can be regarded as 
a direct avenue towards the resolution of these issues - the problem which is currently 
under study. 
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